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Abstract 

We analyze the concept of entanglement for multipartite system with bosonic and fermionic 
constituents and its generalization to systems with arbitrary parastatistics. We use the rep- 
resentation theory of symmetry groups to formulate a unified approach to this problem in 
terms of simple tensors with appropriate symmetry. For an arbitrary parastatistics, we 
define the S-rank generalizing the notion of the Schmidt rank. The S-rank, defined for all 
types of tensors, serves for distinguishing entanglement of pure states. In addition, for Bose 
and Fermi statistics, we construct an analog of the Jamiolkowski isomorphism. 
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1 Introduction 

Confronted with entangled states of a bipartite system, Schrodinger said that entanglement was 
not ...one hut rather the characteristic trait of quantum mechanics, the one that enforces its 
entire departure from classical line of thought. .. [1] . 
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Note that Schrodinger's interests for this aspect of Quantum Mechanics originated from the EPR 
paper. Indeed, the entanglement for states of composite systems created quite an embarrassment 
for theoretical physicists not sharing the point of view of the "Copenhagen interpretation" of 
quantum mechanics. As a matter of fact, the entanglement played a rather important role 
in the development of quantum mechanics, as it obliged the physics community to face the 
nonlocal nature of the description of natural processes. Nowadays, however, entanglement 
is considered as an important resource for quantum computation, quantum information and 
quantum teleportation. 

The usual definition of entangled states is given by defining first separable pure states as de- 
composable tensor products and then declaring a state to be entangled if it is not separable. 

In the case of a system containing identical constituents, the problem of the definition of en- 
tanglement has to be reconsidered, since the indistinguishability produces some amount of 
entanglement by its own. For instance, a non-zero skew-symmetric 2-tensor is never decompos- 
able. We cannot therefore simply transpose to the Hilbert space of the composite system the 
requirement that non-entanglement is guaranteed by the factorization of the total state into 
a simple tensor product of vectors corresponding to the two subsystems. We have to analyze 
better the meaning of entanglement per se and refine the concept of separability. 

Several nonequivalent ways of identification and quantification of the phenomenon were pro- 
posed. In [2j the authors explicitly stated the problem of distinguishing non-local correlations 
from those caused by the Pauli principle measured in fundamental experiments, aiming at 
checking the validity of quantum mechanics involving two identical fermions. To quantify the 
amount of extra correlations, the authors used the number of terms in the decomposition of the 
wave function in terms of elementary 2x2 'Slater determinants' (see Section [6] below) . Similar 
ideas, based on the lengths of the canonical decomposition of the wave function, were used in 
[3] for description of correlations in double ionization of atoms in strong electromagnetic fields. 

A scheme of extracting the correlations caused by the Pauli principle in fermionic systems in 
order to characterize a genuine entanglement, became important in analysis of elementary quan- 
tum gates based on quantum dots [1] . Subsystems here are no longer separated by macroscopic 
distances; on the contrary, they occupy the same spatial regions and their indistinguishability 
becomes relevant. The problem was analyzed in [5] and further in [6j in terms of modern theory 
of entanglement. Measures of correlations were constructed in analogy with the distinguishable 
particles case, again by employing algebraic properties of the coefficient matrix in the expansion 
of a state in terms of basis states. 

Similar ideas can be applied, and in fact were applied, to bosons [H [7]. In the bosonic case 
non-entangled states are again identified with simple symmetric tensors (see Section U] below). 
Thus, in the case of two bosons a state is non-entangled if and only if it is a tensor product of two 
identical one-particle states. This definition should be contrasted with another one advocated 
already in [S], where bosonic and fermionic systems were treated in parallel. To identify two 
partite systems with minimal identical-particle correlations, i.e., without entanglement, one 
considers a situation when two particles are confined to two non-overlapping spatial regions in 
which they are subjected to spatially confined independent measurements. The correlations are 
minimal if, for every independent choice of measurements, the probability of outcomes factorizes 
into the product of single particles probabilities of outcomes. For two fermions, the result 
reduces to the demand that the wave function of the whole system is a simple antisymmetric 
tensor (a single Slater determinant), \^) ~ (8> \4>) — |</>) i-c, to the very definition 
adopted by other authors cited above. For bosons, however, the definition leads to two different 
classes of minimally correlated state. In addition to tensor products of identical one-particle 
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states, also wave functions in the form of Slater permanents i.e., symmetrized tensor products 
of orthogonal one-particle states, I*?) ~ {\ip) (8) |0) + ^ |^)), {tp\(p) = 0, are classified as 
minimally entangled^ The definition can be easily extended to many particle systems [9]. 

In a series of papers by Ghirardi et. al. [lOl [TTl [12], a definition of non-entangled state was 
based on the concept of a complete set of properties possessed by a constituent of a composite 
system. For a bi-partite system, in a pure state p, we say that a one of the constituents has a 
complete set of properties if and only if there exists a rank-one projection operator P acting in 
a single particle space such that Ti{Ep) = 1, where E = P(^I + I(^P — P^P is a projection 
operator acting in T-L A 7i oi 7i V Ti for, respectively, fermions and bosons. This idea resembles 
the one advocated in the above cited papers of Herbut, but rather than invoking directly 'local' 
measurements it stresses properties of subsystem states. The result is the same. For fermions, 
non-entangled states are antisymmetrized tensor products, whereas for bosons they split into 
two classes: products of two identical states, or symmetrized products of two orthogonal states. 
Whereas for fermions the criterion of non-entanglement, based on the number of coefficients in 
expansion in terms of elementary 2x2 Slater determinants, remains functional, it is not the 
case for bosons. Here, states expressible as linear combinations of two Slater permanents are 
non-entangled if the both coefficients of expansion are the same - such states can be expressed 
as symmetrized products of two orthogonal vectors [9l [10]. A concept of entanglement for 
indistinguishable particles based on a measurability of correlations was recently considered also 
in [13]. 

In this paper we attempt to approach this problem from a mathematical view point, where 
structures and available mathematical constructions are used as a guide. A fundamental char- 
acter of our work depends also on reviewing and generalizing basic concepts of mathematical 
foundations in understanding entanglement. We shall provide few physical considerations in 
the conclusions, where we also briefly compare our results with other outlined above. Here, 
let us mention only that our natural and unifying mathematical model strongly suggests non- 
entanglement for bosons to be associated with tensor products of identical states. 

The paper is organized in the following way. In Section [2] we recall definitions and basic 
facts from tensor algebra underlying our analysis of entanglement for multipartite bosonic and 
fermionic systems. In Section [3] we introduce the concept of duality and analyze contractions 
between dual spaces of tensors, which allows us to define the S-rank of a tensor, generalizing 
the Schmidt rank, and then the simplicity of a tensor in Section [H Section [5] contains various 
characterizations of simplicity of a tensor in general, and in the bosonic or fermionic class. In 
Section [6] we define entanglement for bosonic and fermionic multipartite states and provide a 
simple characterization of entanglement for pure states in terms of bilinear functions in coeffi- 
cients of representing tensors. In analysis of entanglement of distinguishable particles, the so 
called Jamiolkowski isomorphism played the role of a very useful tool. We give an extension 
of it for boson and fermions in Section [71 The mathematically rigorous generalization of en- 
tanglement to multipartite system with arbitrary parastatistics, another novel invention in this 
paper, is given in Section [Sj 

Note that such an unifying approach allowed also for description of pure non-entangled states as 
the images of generalized Segre maps, in full analogy with the case of distinguishable particles 
[lll[15]. These questions, however, we decided to postpone to a separate paper. 

^To make the non-entanglement condition in bosonic and fermionic systems completely parallel, one can 
demand orthogonality also in the latter case, since it does not play a role for fermions. 
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2 Tensor algebras 



To describe some properties of systems composed of indistinguishable particles and to fix the 
notation, let us start with introducing corresponding tensor algebras associated with a Hilbert 
space H. For simplicity, we assume that T-L is finite-, say, n-dimensional, but a major part of our 
work remains valid also for Hilbert spaces of infinite dimensions. Note only that in the infinite 
dimensions the corresponding tensor product Hi 'H2 is the tensor product in the category of 
Hilbert spaces, i.e., corresponding to Hilbert-Schmidt norm. 

In the tensor power T-L®^ = ■ - ^H , we distinguish the subspaces: H^*^ = Ti V ■ y V Ti of 

fc-times fc-times 

totally symmetric tensors, and H^'' = Ji f\ ■ - f\'H, of totally antisymmetric ones, together with 

/;■— times 

the symmetrization, tTj^ : "H®*^ H^'', and antisymmetrization, : T-L®'^ T-L^^, projectors: 

T^kUl®---® fk) = ^---'S'/aCfc), (1) 

^^{fi^---®fk) = ;^ E(-l)"/-(i)<^---®/-(fe)- (2) 

■ aeSk 

Here, is the group of all permutations a : {1, . . . , A:} — > {!,..., k}, and (—1)'^ denotes the 
sign of the permutation a. Note that with every permutation a E Sk there is a canonically 
associated unitary operator Ua on Ji^^ defined by 

Uaifl (8) • • • (8) /fe) = fa{l) (8) • • • (8) fa{k) , 

so that the map a 1-^ is an injectivc unitary representation of Sj. UiH®^). We will write 
simply (T instead of U^, if no misunderstanding is possible. Symmetric and skew-symmetric 
tensors are characterized in terms of this unitary action by a{v) = v and a{v) = {—lYv, 
respectively, for all a e Sk 

We put, by convention, = = = C. It is well known that the obvious structure of a 

00 , 

unital graded associative algebra on the graded space = (8 H'^ (the tensor algebra) induces 

k=0 

00 

canonical unital graded associative algebra structures on the spaces "H"^ = © Ti"^ (called the 

k=0 

00 , 

bosonic Fock space) and % = Q) T-L (called the fermionic Fock space) of symmetric and 

fc=o 

antisymmetric tensors, with the multiplications 

V i;2 = 7r^(vi (8 ^^2), (3) 

Wl AW2 = TT^{wi (8 W2)- (4) 

Here, of course. 



= e vr^ : H« ^ H^, (5) 

fc=0 



and 



TT^ = © TT^ : ^ H^, (6) 

k=0 

are the symmetrization and antisymmetrization projections. Note that the multiplication in 
T-L^ is commutative, vi\/V2 = V2V vi, and the multiplication in T-L^ is graded commutative, 
wiAw2 = (-l)^i''^2y;2 A Wl, for Wi G H^'^'. 
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Denote with H* the complex dual space of H. The Hermitian product (•!•) on H induces a 
canonical bijection between H and H* which, in the Dirac's notation, reads 



'H3\x)^ {x\ G n*. 



We must stress that this is not an isomorphism of complex linear spaces, since the above map is 
anti-linear. Note, however, that the symmetric tensor algebra 'Hy can be canonically identified 
with the algebra Pol{T-L*) of polynomial functions on the complex dual %* of %. Indeed, any 
/ G H can be identified with a linear function Xf on H* by means of the canonical pairing 
{,) : Ti X %* —7- C between the dual spaces, as Xf{y) = (/, y) (we must distinguish this pairing 
from the Hermitian product on ^). This can be extended to an isomorphism of commutative 
algebras in which /i V • • • V corresponds to the homogenous polynomial Xf^ ■ ■ - Xf^. Similarly, 
one identifies with the Grassmann algebra Grass{T-L*) of polynomial (super)functions on 
%* . Here, however, with f € Ti we associate a linear function on Ti* regarded as and odd 
function: ^/^/' = — C/'C/- the language of supergeometry one speaks about the purely odd 
manifold HH* obtained from the standard (purely even) linear manifold Ti* by changing the 
parity of linear functions. In this sense, is the algebra of holomorphic (super)functions on 
the complex supermanifold H*. 

If we fix a basis ei, . . . , e„ in ^ and associate with its elements even linear functions xi, • • • , Xn 
on H*, and odd linear functions Ci, • • • , on IIT-L*, then yy ~ C[xi, . . . , i.e. becomes 
isomorphic with the algebra of complex polynomials in n commuting variables. Similarly, ^ 
. . . ,^n]; i-e., is isomorphic with the algebra of complex Grassmann polynomials in n 
anticommuting variables. The subspaces T-L"^^ and Ti^^ correspond to homogenous polynomials 
of degree k. It is straightforward that homogeneous polynomials x'l^ ■ ■ ■ x^", with ki + ■ ■ ■ -\- 
kn = k, form a basis of 7^^^, while homogeneous Grassmann polynomials Cii ^ " " ' Cife) with 
1 < ii < i2 < • • • < ^fc < n, form a basis of H^^. In consequence, dim'H'^^ = {"'^^~^) and 
diuiT-L^^ = (^), so the gradation in the fcrmionic Fock space is finite (for finite-dimensional T-L). 
Of course, we can put together both algebras and consider the tensor product Tii ' with 
dimiJi = n and dimi?2 = rn, which is a graded associative algebra with a bi-gradation N x N 
concentrated in N x {0, 1, . . . , m}. 



Here, T-L^ represent systems composed of k bosons, described by a Hilbert space of 
dimension n, and / fermions in the Hilbert space of dimension m. The whole algebra Tii (8) 7^2 
is the algebra of polynomial functions on the linear supermanifold HI x 11^2 of dimension 
(n,m). Such functions are written as finite complex linear combinations 



nt'^n^'. 



(7) 



ik,l)eNx{0,l,...,m} 





(8) 



ki .^kn 
l<ii< - ■ ■ <il<m 



Note that any basis {ei, . . . , e„} in H induces a basis {ei^ (^Cij®- • •'^ei^, | ii, . . . , ifc G {1, . . . 
in T-L®^. Therefore, any u G l^P^ can be uniquely written as a linear combination 



n}} 



n 




(9) 



«i,...,ifc=l 
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If u G T-L^^ ^ then the tensor coefficients n*!---** are totally symmetric and, after applying the 
symmetrization projection to Q, we get 

n 

u= u'^-^'e^.V ...Vei^. (10) 

il,...,ik = l 

Similarly, if u G Ti^^, the tensor coefficients n^^" **^ are totally antisymmetric and 

n 

u= u''-'>'ei,A...Aei^. (11) 

ii,...,ifc=l 

We will refer to the coefficients u*^"'** as to the coefficients of u in the basis {ei, . . . ,e„}. 



3 Tensor duality and contractions 

Starting with the canonical duality (pairing) ( , ) between vectors from Ti and covectors from 
H*, with an obvious prolongation to a paring between "H®^ and {T-L*)®^ , 

k 

= J\{fhgi), (12) 

i=l 

and viewing symmetric and antisymmetric tensors as canonically embedded in the tensor alge- 
bra, we find canonical pairings between T-L^^ and {'H*)'^^, as well as between T-L^^ and {T-L*)^^ . 
For /i, . . . , /fc G and gi, . . . , gk G V.* , we get 

k 

(/iV--- V/fe,<7iV---Vgfc) = — -2 n(/-(i)'9r(i)> = ;^per((/i,5i)). (13) 

Here, X^^g^^, 0^=1 ^iT(j) = per(ajj) is the permanent of the matrix A = (aij). Similarly, 

(/i A--- A/fc,5i A--- Agfc) = ^det((/„(7j)). (14) 

Quite similarly one can prove that the natural Hermitian product on T-L^, for which the grading 
is the decomposition into orthogonal subspaces and 

k 

{fl(^---(^fk\f[^---C^f'k)= ll(fi\f^)^ (15) 

1=1 

induces Hermitian products on the subspaces H"^^ and which read, respectively, 

(/iV---VM/{V---V/^) = lper((/,|/j)), (16) 
(/iA---AM/l A---A/^) = ldet((/,|/j)). (17) 

Given a basis fi, . . . , fn of Ti and the dual basis f^, . . . , f* of Ti* , we have the induced bases: 

f^' V • • • V , ki + --- + kn = k 
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of and 

fii ^ fik^ 1 < n < Z2 < • • • < ■jjt < n 
of H^^. The dual bases read 

h\.'.K\ ' fcl + --- + ^n = fc, (18) 

and 

fcl/fi A---A/*, l<^l<^2<•••<^ik<n. (19) 
Analogously, any orthonormal basis ei, . . . , e„ of induces canonical orthonormal bases 



^-y-^^ V---Ve^, ki + --- + kn = k, (20) 
of -^v^, and 

Vklei^ A • • • A Cj^, 1 < < i2 < • • • < ifc < n , (21) 

of n^''. 

The canonical pairings between 'H^'' and {%*Y^ on one hand, and and (^*)'^'^ on the 
other, can be generalized to certain 'pairings' [contractions or inner products) between T-C'^ 
and {T-L*Y^ on one hand, and H^*^ and {H*)^^ on the other. For the standard simple tensors 
/ = /i <8 • • • <8 /it e H®^ and 5 = 51 ® • • • ® 5i G (H*)®^ we just put 

^gf = {fi® ■■■ ^ fl,9i^ ■■■ ® 9l)fl+i ^■■■^ fk 

if I < k and igf = if / > /c, and extend it by linearity to all tensors. It is easy to see now that, 
if v = fiV ---V fkeH"^'' C and = 51 V • ■ • V 5, G {H*Y^ C (H*)®' then i^v G H"^^^-^\ 

Similarly, i^w G 'H^^^-^\ if u; G H^*^ C U®^ and a; G (H*)^' C (H*)®'. Explicitly, 

1 ' 

' o-eSk j=i 

reSi 

= ^-^jr^ E U(fsu)^srij))fsm) V • • • V fsik), (22) 

res, 

where S{l,k — I) denotes the set of all {l,k — I) shuffles. 
For skew-symmetric tensors, 

1 ' 

IgiA-AgJl ^ ' ' ' ^ fk = i''^)" (''^T Yl{fa{j), 9T(,j)) fa{l+l) <8 " " " ® fa{k) 

' TeSk j=i 

= ^V^ E (-ir(-irri(/5o>5.o-))/5(w)A---A/5(fe). (23) 

SeS(l,k-l) j=l 
reSi 

In particular, 

W-Vfl Ji V • • • V /fe = (/i V • • • V /fc, 51 V • • • V 5fe) = ^per((/i |5,)) , (24) 
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and ^ 

hiA-AgJi A • • • A /fc = (/i A • • • A A • • • A = — det{{fi\gj)) . (25) 

Moreover, 

-^^ k j 

Vv-vp,_i/iV---V/fe = V---V5fe_i)/j (26) 

1 

and 

-[^ i 

%iA-A3,_i/i A--- A/fc = ^^(-l)^'"^(/i A--- A/fc,5i A--- A(7fc_i)/,- (27) 

fc 



i^(-l)^-^det((/,|5.).^,)/,. 

" i=i 



4 The S-rank 



There are many concepts of a rank of a tensor used to describe its complexity. One of the 
simplest and most natural is the one based on the inner product operators defined in the 
previous section. We will call it the S-rank, since it turns out to be a natural generalization of 
the Schmidt rank of 2-tensors. 

Definition 4.1. Let u G U®^ . By the S-rank of u, we understand the maximum of dimensions 
of the linear spaces i^^a{u), for a E S^, which are the images of the contraction maps 

3 ^ ^ g ^ (28) 

Theorem 4.1. If u e T-C^ (resp., u G T-L'^^), then the S-rank of u equals the dimension of the 
linear space being the image of the contraction map 

(^*)V(fe-i) 3 ^ ^ g (29) 

(resp., 

{n*)^^^-"-^ 3 li^i^ueU). (30) 

Proof: It follows immediately from the observation that a contraction of a symmetric (res., 
antisymmetric) tensor u with a tensor ^ G (H*)®' is the same as its contraction with the 
symmetrization (resp., antisymmetrization) of /x, z^n = It^\/(^^^u, (resp. i^u = «7r^(/i)'") that, 
for a G Sk, a{u) = ±u. 

□ 

Example 4.1. Let ei,e2 be orthogonal vectors in 7i. The vector 

n = ei V 62 = ^ (ei (g) e2 + 62 (8) ei) (31) 
has S-rank 2. Indeed, for each x eH, 

ixU = (ei 62 + 62 (8) ei) = ^ ((x|6l)62 + (x|62)6l) , 

so that ly^u is spanned by {61,62}. 
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Example 4.2. Let us take T-L of dimension 3 and an orthonormal basis Cj, i = 1,2,3. Let 
u G T^y^ in the polynomial notation reads 

u = + e^ + 16e^ + 4e?e2 + 8e?e3 + 4eiei + 86^63 + 32eiei + 32e2ei (32) 
+66^62 + 24eie| + 2Aelel + 24616263 + 24616263 + 486i626| . 



The vectors 



3332 2 22 2 2 

ei > 62 , 63 , 6162 , 6163 , 6162 , 6263 , 6163 , 6263 , 616263 



form a basis of T-C^^ and it can be directly checked with the use of (|26|) that any of these vectors, 
say /U, contracted with u gives a vector proportional to 61 + 62 + 263. For instance, the only 
non-zero parts of i^u for ^ = are 

i^i{u) = Zg3 (62 + 46162 + 86263) = 62 + 61 + 263 

and 

lelesiu) = Z,2,3 (86? 63 + 246? 61 + 246^6263) = ||e?63f (661 + 1263 + 662) = 2(6i + 263 + 62) ■ 

It follows that iy^vs{u) is 1-dimensional, so the S-rank of n is 1. 

Example 4.3. For Ti with the basis as above, consider w G Ti^^ of the form 

t(; = 61 A 62 + 62 A 63 . (33) 
Let us see that the S-rank of w is 2. Indeed, 



1 ^ 1 

-62, ie2W = - [es - ei) , i^^w = --( 



so that i-uw is spanned by 62 and 63 — 61. 

Example 4.4. For T-L with an orthonormal basis 61,62,63,64 the tensor 

w = 61 A 62 + 63 A 64 (34) 

has the S-rank 4. Indeed, 

1 11 1 

^ei""^ = 2^2 , ie2W = --ei, le^W-Ci, ieiW = --63 . 

Theorem 4.2. (a) The minimal possible S-rank of a non-zero tensor u G equals 1. A 

tensor u G Ji®^ is of S-rank 1 if and only if u is decomposable, i.e., u can be written in 
the form 

u = fi®---®fk. fi^n, fi^O. (35) 

Such tensors span T-L®^ . 

(b) The minimal possible S-rank of a non-zero tensor v G 'hy^ equals 1. A tensor v G 'hy^ is 
of S-rank 1 if and only if v can be written in the form 

v = fv---vf, fen, //O. (36) 

Such tensors span . 
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(c) The minimal possible S-rank of a non-zero tensor w S Ti. equals k. A tensor w G T-L is 
of S-rank k if and only if w can be written in the form 

w = fiA---Afk, (37) 

where fi,...,fk € are linearly independent. Such tensors span H^^. 

Proof: The 'if parts of the above statements are obvious, so we shall prove 'only if. If the 
dimension of t'^^u is 1, thus the space is spanned by some f = fk £ Ti, then clearly u = u' ^ f 
for some u' G T^^C^-i). If^ in turn, the tensor is symmetric, then clearly u is proportional to 
/(g) • • • (g) /. In the general case we get a similar fact for a{u) with a ^ Sk, so u = fi® ■■■ ® fk- If 
u is skew-symmetric and fi, ■ ■ ■ , fr span ly^^u, then w is a linear combination of tensor products 
of these vectors, so of /j^ A • • • A /j^, . Hence, r > k, and r = /c if and only if u is proportional 
to /i A • • • A /fc. That simple tensors span the corresponding spaces is pretty well known for 
general and antisymmetric tensors. For symmetric tensors it follows from the fact that powers 
of linear functions span the spaces of polynomials with coefficients in a field of characteristic 0. 

□ 

Definition 4.2. Tensors of minimal S-rank in Ji^^ (resp., 7^^'^, Ti^^) we will call simple (resp., 
simple symmetric, simple antisymmetric). 

Theorem 14.21 immediatelv implies the following. 

Corollary 4.1. The S-rank is 1 for simple and simple symmetric tensors, and it is k for simple 
antisymmetric tensors from TL^^. Simple tensor have the form i35]) . simple symmetric tensors 
have the form i36\) . and simple antisymmetric tensor have the form |g7| j. 

Example 4.5. The symmetric tensor u defined in (|32p . Example 14.21 has the S-rank 1, so it is 
simple. As a matter of fact, 

u = {ei + e2 + 263)'' . 
Also the antisymmetric tensor w defined in (j33|) . Example 14. 3[ is simple: 

= 62 A (63 - ei) . 

Remark 4.1. Of course, for distinguishable particles there is no need to use the same Hilbert 
space 7i in the tensor products. We can use the tensor product Tii g) • • • (g) Tik instead, with 
simple tensors being decomposable: u = /i (g • • • (g /fc, fi G Tii. 



5 Various characterizations of simple tensors 

For n^^, where Uq = Ui ® ■ ■ ■ ® Uk, denote with ai : TL^"^ i = l,...,k, the 

transposition with respect to the i-th. and the {k + i)-th factor, 

ai{ai g) • • • (g a2k) = (38) 
ai (g • • • (g Ui-i g) Qk+i g) aj+i (g • • • g) fflfc+i-i (g aj (g ak+i+i g) • • • (g a2A: , 

and with : T-Lq'^ — )■ T-Lq'^ - the cyclic permutation that moves the last factor into the first 
place: 

rfe(ai g) • • • (g a2k) = a2fc ai g) • • • g) a2k-i ■ 
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Theorem 5.1. For u £ Tii (i^ ■ ■ ■ ^ %k the following are equivalent: 
(a) u is simple; 

(h) Vcj G 5fc V /XI, /i2 G •^^(i) ^ V2f^H = 0; 

(c) ai{u u) = u(d) u for i = 1, . . . ,k. 

Proof: (a) =^ (b) is obvious in view of Theorem 14.21 (a). Also (b) =^ (c) is clear, so assume 
ak{u0u) = ui^u. This implies that the dimension of the space ly^^u is 1. Indeed, if this space 
is spanned by linearly independent vectors gi, . . . ,gr, then u = X^j=i Ur <X' gr for some linearly 
independent uj G T^^C^-i). Since ak{u®u) = u®u means that 

r r 

Uj ® gj ® Us® gs = ^ Uj ® gs®Us® Qj , 

j,s=l j,s=l 

we conclude that r = 1, so u = u' ® fk for some fk^'H and u' E U^i^-i) ^ ^ similar reasoning, 
applied to the identities (Ti{u ®u) = u®u with i = 1, . . . , /c — 1, implies that u = fi® ■ ■ ■ ®) fk- 
□ 

Theorem 5.2. For a symmetric tensor v G T-L^^ the following are equivalent 

(a) V is simple symmetric, 

(b) V z^i, 1/2 G {n*Y^^-^) t^,v A i^^v = 0, 

(c) V ® V = ak{v ® v). 

Proof: (a) =^ (b) is obvious in view of Theorem 14.21 (b). As (b) implies v = v' ® f^ for 
some fk S H, also (b) =^ (c) is clear. The condition (c), in turn, for symmetric tensor yields 
V (gi V = ai{v ® v) ior all i = 1, k, so V = fi ® ■■■ ® fk as above, thus v = f ® ■ ■ ■ ® f hy 
symmetry. 

□ 

Theorem 5.3. For an antisymmetric tensor w G the following are equivalent 

(a) w is simple antisymmetric, 

(b) y cue ujAi^w = 0, 

(c) {T^k+l 'i'dy^S,{k-l)){Tk{w ® W)) = 0. 

Proof: (a) =^ (b) is obvious in view of Theorem l4.2l (c). Also (b) =^ (a) is clear and well known. 
We shall show that (b) and (c) are equivalent. 

Let us write w as a sum of simple antisymmetric tensors, w = X^jLi "W^j; with minimal m. Write 

Wj = fj A ■ ■ ■ A fj and denote Wj = fj A ■ ■ ■ f- A ■ ■ ■ A fj , s = 1, . . . ,k. Here, stands for the 
omission. As the number of simple tensors is minimal, the tensors Wj are linearly independent 

in w^^-^y 

Observe now that 

s 
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so that 

\k—s 



( — 1 ) 



Since Wj are linearly independent, the latter vanishes if and only if all fj A w vanish, that is 
clearly equivalent to (b). 

□ 

Remark 5.1. The conditions (c) in Theorems 15.11 [5^21 and 15.31 have formally this advantage 
over the corresponding conditions (b) that they are directly verifiable, as they do not contain 
general quantifiers referring to infinite sets. 

Note that all we have said remains valid for an arbitrary vector space over a field of character- 
istics 0. The Hermitian structure played no role yet. 

Let us note first that the S-rank is associated with a point in a projective space rather than with 
a tensor itself. Hence, we can restrict considerations to tensors of length 1 (as vectors in the 
Hilbert space Hi ® ■ ■ ■ (^lik)- According to Theorem 15.11 a tensor u £ Tii ■ ■ ■ Tik simple 
if and only if u0 u = ai{u<^ u) for all i = 1, . . . , A:, where <Tj interweaves the two copies of Hi in 
^1 (g) • • • (g) T-Lk (g) (g) • • • (g) Hk- Note that fij acts as a unitary operator in {Hi (g • • • (g) H.k)'^'^ ■ 
This means, in turn, that 

k 

^\\u®u-ai{u®u)\\^ = ^. (39) 
1=1 

The latter we can write as 

k 



(g) u — ai{u ®u)\u®u — ai{u ® u)) = 0, (40) 

1=1 

which, for tensors of length 1 is equivalent to X^(Li ^^{u (g) u\(Ti{u (g) u)) = k, or, finally, to 

Re (g) u I (n (g) u)^ = k . 

The Schwarz inequality yields now u0 u = a{u (g u), where o' = ^ Y17=i 



1 

a{xi 1^ ■ ■ ■ (S> Xki^ yi 1^ ■ ■ ■ 1^ Vk) = - y^^xi(S> ■ ■ ■ 1^ ■jji'^ ■ ■ ■ yi'S> ■ ■ ■ Xi(Si ■ ■ ■ (gi yk ■ 

i=l 

In this way we have proven the following. 

Theorem 5.4. Let u G Hi (g) • • • (g) Hk, \\u\ \ = 1- Then, u is simple if and only if 

{u (g u, a{u (g) u)) = 1 . 

For symmetric tensors, a similar fact can be proven analogously. 
Theorem 5.5. Let u G , ||n|| = 1. Then, u is simple symmetric if and only if 

{u g) n, ak{u (g) u)) = 1 . 
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The structure of a Hilbert space % has been used in [H El [6] to define simple symmetric 
and simple antisymmetric tensors for tensors of order 2 by means of the Slater decomposition. 
Both decompositions are direct consequences of a possibility to write a complex symmetric 
(antisymmetric) matrix in a diagonal (block-diagonal) form by a unitary change of basis (Takagi 
theorem [16j). The existence of the Slater decompositions means that any symmetric tensor 
V €z Hy T-L and any antisymmetric tensor w ^ H AH can be written as 

n 

V = ^AjBiVei, Xi>0, (41) 
1=1 

n 

w = ^Xifi/\fn+i, Ai>0, (42) 
1=1 

for some orthonormal systems (e^) and (/j) of T-L. The Slater rank is the number n of terms 
in these decompositions. The above are clearly symmetric and antisymmetric analogues of the 
Schmidt decomposition: any tensor u G Tii 'H2 can be written in the form 

n 

ti = J^AiBiO/i, Ai>0, (43) 
1=1 

for some orthonormal systems: (e^) in Tii, and (/j) in T-L2- 

Theorem 5.6. For any 2-tensor u E Hi (X" ^2 its Schmidt rank equals its S-rank. 

Proof: It is clear that the Schmidt decomposition ()43p implies that the S-rank of u is n. 
Conversely, if the S-rank of u is n, then u can be written in the form (|43p with the only 
difference that the systems (cj) and (fi) are merely linearly independent. But this is a standard 
procedure, used in the proof of existence of the Schmidt decomposition, that we can choose (e^) 
and (fi) orthonormal. 

□ 

Hence, the 2-tensors are simple (resp., simple symmetric, simple antisymmetric), if there exists 
a Schmidt (resp. Slater) decomposition with a single Aj > 0, i.e. they have the Schmidt (Slater) 
rank 1. Unfortunately, there are no direct analogues of these decompositions for tensors of 
higher orders. 

On the other hand, as we have already seen, such type of a decomposition is not necessary to 
define (and check) which tensors are simple symmetric (resp., simple antisymmetric), as the 
S-rank can serve in these cases. 



6 Entanglement for multipartite Bose and Fermi systems 

Using the concept of simple tensors we can define simple (non-entangled or separable) and 
entangled pure states for multipartite systems of bosons and fermions. 

Definition 6.1. 

(a) A pure state p^, on Ti'^'' (resp., on V.^^), = j^^, with v G U^^ (resp., v £ U^^), 
w 7^ 0, is called a hosonic (resp., fermionic) simple (or non- entangled) pure state, if v is 
a simple symmetric (resp., antisymmetric) tensor. If v is not simple symmetric (resp., 
antisymmetric), we call a bosonic (resp., fermionic) entangled state. 
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(b) A mixed state p on 'hy (resp., on T-L ) we call bosonic (resp., fermionic) simple (or 
non- entangled) mixed state, if it can be written as a convex combination of bosonic (resp., 
fermionic) simple pure states. In the other case p is called bosonic (resp., fermionic) 
entangled mixed state. 

According to Theorem 14.21 bosonic simple pure A:-states are of the form 

|eV--- Ve)(eV-- - V e| 
for some unit vector e E "H, and fermionic simple pure /c-states are of the form 

k\\ei A • • • A efc)(ei A • • • A e^l 
for some orthonormal system ei, . . . , in T-L. 

Fixing a base in Ti results in defining coefficients •••**:] of n G Ti^^. Formulae characterizing 
simple tensors, thus simple pure states, can be written in forms of quadratic equations with 
respect to these coefficients as follows. The corresponding characterization of entangled pure 
states are obtained by negation of the latter. 

Theorem 6.1. The pure state pu, associated with a tensor u = G T-L'^^ , is entangled if 

and only if there exist ii, . . . ,ik, ji, ■ ■ ■ , jk, and s = 1, . . . , k such that 

Proof: The tensor u(^u has coefficients u^'^—^ky^ji—kk^ gg Eg . (|44p expresses the fact that u(^u 7^ 

CTg {u®u), and thus Theorem 16.11 is a direct consequence of Theorem 15.11 

□ 

Theorem 6.2. The bosonic pure state p^, associated with a symmetric tensor v = £ 
is bosonic entangled if and only if there exist ii, . . . ,ik, ji, ■ ■ ■ , jk, such that 

yil---ik-iikyjl---jk-ljk yh---ik-ijkyji---jk-iik ^45^ 

Proof: A direct consequence of Theorem 15.21 
□ 

Example 6.1. The bosonic pure state associated with the symmetric tensor u defined in (|31|) 
is entangled, as we know it has the S-rank 2. To apply our characterization ()45p . note that the 
only non-zero coefficients are u^^ = n^^ = ^. Hence, 

n 11 22 / 12 21 ^ 
= u u u u =-. 

4 

Example 6.2. The bosonic pure state associated with the symmetric tensor u defined in (j32|) is 
non-entangled, as we know it has the S-rank 1. To apply our characterization ()45p . note that, as 
easily checked, u*^*^*^**' = 2*(^\ where i(3) is the number of 3's in the sequence i = (^i, ^2) ^3; u)- 
Now it is clear that 

yili2i3iiy^jij2j3ji — 2*(3)+i(3) 

where i{3) +i(3) is the number of 3's in the sequence (ii, ^2, ^a, Hi ji, J2i JSi j4) which remains 
unchanged under permutations. In particular, 
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Theorem 6.3. The fermionic pure state pw, associated with an antisymmetric tensor w = 
Jy;«i---«fc] g Ji^^ J is fermionic entangled if and only if there exist ii,. ■ ■ jifc+i, ji, • • • ,jk-i such 
that 

^[h-ik^ik+i]h-jk-i -^0, (46) 

where the left-hand side is the antisymmetrization of w^^'"^'=w^''+^^'^'"^'=-^ with respect to indices 
ill ■ ■ ■ h+i- 

Proof: A direct consequence of Theorem 15.31 
□ 

In view of the above characterizations, it is obvious that the sets of entangled (entangled bosonic, 
entangled fermionic) pure states are open: pure states sufficiently close to entangled ones are 
entangled. 

Example 6.3. The fermionic pure state associated with the antisymmetric tensor w defined in 
(|33p is non-entangled, as we know its S-rank is 2. To apply our characterization ()46p . let us note 
that the only non-zero coefficients are w^"^ = —w"^^ = w"^^ = —w^"^ = ^. The antisymmetrization 
y^[n«2-y,«3]i (jg^]2 uou-zcro only if (ii,i2)^3) is a permutation of (1,2,3). But, as easily seen, 

^[12^3]i ^ 1 r^u 3j _ 32 lj^ ^ 
3 ^ ^ 

for j = 1,2,3. 

Example 6.4. The fermionic pure state associated with the antisymmetric tensor w defined 
in (|34p is entangled, as we know its S-rank is 4. To apply our characterization ()46p . let us note 
that the only non-zero coefficients are w^"^ = —w"^^ = w'^^ = —w'^^ = ^. As easily seen, 

7 Jamiolkowski isomorphisms for bi-partite systems of bosons 
and fermions 

In the theory of entanglement there exists a useful tool for investigating, on one hand, entangle- 
ment properties of states and, on the other, structure of positivity-preserving and hermiticity- 
preserving maps between matrix algebras called the Jamiolkowski isomorphism |T7]. Origi- 
nally proposed in [TS] as an instrument for checking the property of preserving positive semi- 
definiteness for a linear map between two matrix algebras in finite dimensional spaces, it was 
later used to prove the so-called operator form representation for linear maps on quantum 
states (i.e. positive semi-definite operators) [19]. The construction can be extended to the 
infinite-dimensional setting if the restriction to Hilbert-Schmidt operators is imposed [20j) 
which significantly heightens the usefulness of the Jamiolkowski isomorphism in cases when 
finite-dimensional description of q quantum system is not suitable. An in-depth description 
of the applicability the Jamiolkowski isomorphism is given in |21j where various properties of 
entanglement and separability of states are paralleled with features of the corresponding linear 
maps. 

One of the most straightforward applications of the Jamiolkowski isomorphism is a characteri- 
zation of separability of a state in terms of the rank of the corresponding map. In the following 
we show how this decription can be extended to states of indistinguishable particles. 
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Let us recall (see e.g. |20j) that for two Hilbert spaces Tii and we have the following diagram 
consisting of Jamiolkowski isomorphisms: 



J 



■ C2{C2{'H2,ni)) 



(47) 



Jl 



J2 



£2(^1 ® ^2) 

Here, with C2{'Hi,7i2) we denote the Hilbert space of Hilbert-Schmidt maps from Tii into 'H2- 
Of course, this space reduces to all complex linear maps, if Tii or ^2 is finite-dimensional. Note 
that we write shorter C2{T-L) for £2(^1 Note also that according to the obvious identification 
£2(^2, "Hi) ~ "Hi (8) the identification 



c{'H2,'Hi)* ^c{ni,n2) , 

induced by the natural pairing 

£ (%, ni) X £ (7^1, 3 {A, B) ^ Tr{A o B) e C , 
we can rewrite the above diagram in the form 

J 



ni(^ni®n2(^n*2 



'Hi®n*2(^n2®ni 



(48) 



(49) 



(50) 



Jl 



J2 



ni®n2®n*2(^'H\ 

in which the Jamiolkowski isomorphisms reduce to appropriate permutations of tensors. In the 
case when Hi = ^2 the whole picture reduces to 



t^2{C2{n)) 



J 



■£2(£2(-H)) 



(51) 



Jl 



Jl 



ii2{n®n) 

We can now decompose 11®% into {V. V V.) ® {Ti A V.) and reduce our diagram to maps on 
each of these components. In this way we get 




(52) 



£2(-H^')e£2(^^') 
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Here, Cf-{Ti®'H.*) is the space of self-adjoint Hilbert-Schmidt operators <I>, $ = on H^'H* = 
C2{%) according to the identification C2{%) — {C2{'H))* (see (|^ ) related with the pairing (jMj) . 
Note that this is different from hermicity, <I>* 7^ <I>^, since depends on <^ anti-linearly. 

Indeed, $ written in the Dirac notation as <I> = \ijki\^i) ® \ ® \&k) ® is self-adjoint if 
and only if Tr(A o ^{B)) = Tr($(^) o B) for ah A,B ^ C2{'H), that applied to A = \ej) (g) {ei\ 
and B = \ei) 1^ {ckl yields \ijki = XkUj-, so that we deal with maps coming from symmetric or 
antisymmetric tensors. What is more, since 

J (|ei) (g) {ej\ (g) \ek) {ei\) = \ei) (g (e^l (g \ek) (g {ej\ , 

we have a further splitting 

cfin n*) = £f +(?^ ® ?^*) © £f "(7^ (g n*) , 

where Cf^^iT-L^T-L*) consists of these ^ for which J{^) = ±$. Finally, we end up with bosonic 
and fermionic Jamioikowski maps: 



ci^+in^n*) ^ci^+{n®n*) (53) 



and 




If now p = \v){v\ is a pure state in "hy^ corresponding to a vector v G "H^^ with a Slater 
decomposition v = '^^^^ ^ -^^ ^ 0' ^'^^^ Slater rank is r, then p = J']^(<I>), with 

<I> = ^ \i\j\ei) (g (ejl (gi |ei) (g (ej| 
being a map from C^^iH (g) ?^*) of rank r"^ . 

Similarly, li p = \w){w\ is a pure state in "^^^^ corresponding to a vector G "^^"^^ with a Slater 
decomposition w = Yl\=i f^ifi ^ fn+i, fJ-i > 0, so that the Slater rank is r, then p = J'^{^), with 

$ = ^ PiPj ilfi) © {fj\ © «) {fn+j\ - \fn+i) «> (/il ©!/»)«> (/n+jl 

- I/,) © (/n+jl «) l/n+i) © (/jl - \fn+i) © (/n+j| © «> 

bemg a map from C^'^'i'H © ?^*) of rank 4r2. In this way we get a characterization of bosonic 
and fermionic simple pure states in terms of the corresponding Jamioikowski isomorphisms. 
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Theorem 7.1. A pure state p in'H^'^ (resp., T-L^'^) is bosonic (resp., fermionic) simple if and 
only if p = J+{<^) (resp., p = Jf{^)) for $ e £f +(?^ n*) (resp., <^ G £f-(^ n*)) of 
rank 1 (resp., 4)- 

Remark 7.1. Of course, choosing a basis in T-L we can represent the above maps by matrices 
and gets the Jamiolkowski isomorphism in the form of permutation of indices, more famihar 
to physicists. The above form has the advantage that it does not depend on the basis, i.e. is 
canonical and covariant. 

8 Entangled states of composite systems with generalized paras- 
tatistics 

Our approach to the entanglement of composite systems for identical particles is so general and 
natural that it allows for an immediate implications also for generalized parastatistics. 

Observe first that simple tensors of length 1 in T-L = 7ii® ■ ■ ■ ®%k form an orbit of the group 
U(T-Li) X • • • X U(T-Lk) acting on T-L in the obvious way. In fact, each such tensor can be written as 
for certain choice of orthonormal bases e^, . . . , eii- in Tij, j = 1, . . . , k. This means 
that simple tensors are just vectors of highest (or lowest - depending on the convention) weight 
of the compact Lie group U{(Hi) x • • • x U{T-Lk) relative to some choice of a maximal torus and 
Borel subgroups. If indistinguishable particles are concerned, the symmetric and antisymmetric 
tensors in Ti'^ form particular irreducible parts of the 'diagonal' representation of the compact 
group U{T-L) in the Hilbert space Ti®^, defined by 

U{xi • • • (g) Xfc) = U{xi) (g) • • • U{xk). (55) 

Recall that we identify the symmetry group Sk with the group of certain unitary operators on 
the Hilbert space Ti^ m the obvious way, 

a{xi ® ■ ■ ■ Xk) = (g) • • • (g> x„(k) ■ 

Note that the operators of intertwine the unitary action of U{T-L). In the cases of the 
symmetric and antisymmetric tensors, we speak about Bose and Fermi statistics, respectively. 
But, for k > 2, there are other irreducible parts of the representation (j55p . associated with 
invariant subspaces of the S'fc-action, that we shall call (generalized) parastatistics. Any of these 
/c-parastatistics (i.e. any irreducible subspace of the tensor product Ti'^'') is associated with a 
Young tableau a with /c-boxes (chambers) as follows (see e.g. [221123]). 

Consider partitions of k: /c = Ai + • • • + A^, where Ai > • • • > A^ > 1. To a partition 
A = (Ai, . . . , \r) is associated a Young diagram (sometimes called a Young frame or a Ferrers 
diagram) with A, boxes in the ith row, the rows of boxes lined up on the left. Define a tableau 
on a given Young diagram to be a numbering of the boxes by the integers 1, . . . ,k, and denote 
with Y\ the set of all such Young tableaux. Finally, put Y(k) to be the set of all Young tableaux 
with k boxes. Given a tableau a £ Y{k) define two subgroups in the symmetry group Sk- 

P = Pa = {(7 £ Sk ■ o- preserves each row of a} 

and 

Q = Qa = {c G 5fc : 0" preserves each column of q} . 

In the space of linear operators on Ti'^^ we introduce two operators associated with these 
subgroups: 

a„ = ^a, 6„ = ^ (-1)^(7. (56) 

o-GP ctGQ 
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Finally, we define the Young symmetrizer 



Qq, o be 



E ( 



(57) 



It is well known that vr" = j^Ca, for some non-zero rational number /u(a), is an orthogonal 
projector and that the image Ti" of Cq is an irreducible subrepresentation of U(T-L), i.e. the 
parastatistics associated with a. As a matter of fact, these representations for Young tableaux 
on the same Young diagram are equivalent, so that the constant fj,{a) depends only on the 
Young diagram A of a (does not depend on the enumeration of boxes), /x(a) = /i(A), and is 
related to the multiplicity m(A) of this irreducible representation in Ji^^ by /u(A) • m(A) = kl. 
For a given Young diagram (partition) A, the map 



/.(A)^ 



(58) 



is an orthogonal projection, called the central Young symmetrizer, onto the invariant subspace 
being the sum of all copies of the irreducible representations equivalent to that with a paras- 
tatistics from Yx. 

The symmetrization vr"^ (antisymmetrization vr^) projection corresponds to a Young tableau 
with just one row (one column) and arbitrary enumeration. It is well known that any irreducible 
representation 71°^ of U(T-L) contains cyclic vectors which are of highest weight relative to some 
choice of a maximal torus and Borel subgroups in U{T-L). We will call them a-simple vectors 
or simple vectors in 'H°'. Note that such vectors can be viewed as generalized coherent states 
[23]. They were also regarded as the 'most classical' states by several authors [25]. These are 
exactly the tensors associated with simple (non-entangled) pure states for composite systems 
of particles with (generalized) parastatistics. This is because a-simple tensors represent the 
minimal amount of quantum correlations for tensors in the quantum correlations forced 
directly by the particular parastatistics. 

Example 8.1. (a) For k = 2 we have just the obvious splitting of "H®^ into symmetric and 
antisymmetric tensors: H^'^ ©^^^. 

(b) For A; = 3, besides symmetric and antisymmetric tensors, we have two additional irreducible 
parts associated with the Young tableaux 



respectively. Hence, 
with 



and 



1 


2 


CO 





and Q2 



1 


3 


2 





^03 ^ ^A3 ^ ^ai ^ ^ ^V3 ^ 



(59) 



(60) 



(61) 



-{Xi (g) X2 (S" X3 + X2 (S" Xl (g) X3 - X3 (g) X2 (8) Xl - X3 (g) Xl (g) X2), (62) 
o 



^a2 .-^®3 ^^a2 



(63) 



7r°2(xi (g X2 (g) X3) 



-(xi X2 X3 + X3 (g) X2 (g Xi - X2 Xi (g X3 - X2 X3 (g Xi). (64) 
o 
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The simple tensors (the highest weight vectors) in T-L"^ can be written as 

= A(ei (g> ei (g) 62 - 62 (8) ei (g) 6i), (65) 

for certain choice of an orthonormal basis ei in % and A 7^ 0. Analogously, the simple tensors 
in in turn, take the form 

= A(6i (g) 62 61 - 62 O 61 (g) 61). (66) 

For dim(^) = 3, the simple tensors of length 1 form an orbit of the unitary group U{T-L) of the 
(real) dimension 7 in 7^"^ and Ti"^ . The simple symmetric tensors of length 1 form an orbit of 
the dimension 5, and the simple antisymmetric ones (of length 1) ~ an orbit of the dimension 
1. The dimensions of the irreducible representations are: dim(?^^^) = 1, dim(?^^^) = 10, 
dim(-H"i) = dim(?^"2) = 



Let y." C "H® denotes the irreducible component of the tensor representation of the unitary 
group U{'H) in Ji'^^ associated with a Young diagram a G y{k). 

Definition 8.1. 



(a) We say that a pure state p G 'H^'^ obeys a parastatistics a G Y{k) (is a pure a-state in 

short) if p is represented by a nonzero tensor v G Ti", i.e. 

f = P. = |^^. «i7) 

I Pl 1 

In other words, p is a pure state on the Hilbert space T-L°'. 

(b) A pure state p G Ti^^, obeying a parastatistics a is called a simple pure state for the 
parastatistics a {simple pure a-state, in short), if p is represented by an a-simple tensor 
in Ti'^. If p is not simple a-state, we call it entangled pure a-state. 

(c) A mixed state p on 'H°' we call a simple (mixed) state for the parastatistics a (simple 

a-state in short), if it can be written as a convex combination of simple pure a-states. In 
the other case, p is called entangled mixed a-state. 



For an arbitrary parastatistics (Young tableau) a G Y{k) with the partition (Young diagram) 
A = (Ai, . . . , Xr), consider the map 

ia-'W''^ n^'' , {xi,...,Xr) ^ Xa(l) (g • • • (g X^^k) , 

where a{i) is the number of the raw in which the box with the number i appears in the tableaux 
a. In other words, we make a tensor product of k vectors from {xi, . . . , Xr} by putting xj in the 
places indicated by the number of the boxes in the jth row. For instance, the Young tableaux 
from Example 18. II give iai (xi, X2) = xi (g xi (g X2 and {^1^x2) = xi (g X2 <g xi. It is clear that 
iaixi, . . . , Xr) is an eigenvector of a^. 

It is easy to see now that the S-rank of the tensor 7r"{xa(i) (g • • • (g Xa(k))j for xi A • • • A x^ / 0, is 
r and that this is the minimal S-rank for tensors from "H". Hence, the minimality of the S-rank 
is a good characteristic also for simple a-tensors. 

Theorem 8.1. A tensor u G 7^" is a-simple if and only if u has the minimal S-rank among 
all non-zero tensors from Ti". This minimal S-rank equals r - the number of rows in the 
corresponding Young diagram. 
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9 Conclusions 



States of identical particles exhibit a priori correlations caused merely by (anti)symmetry of the 
wave function in the case of fermions or bosons. It is thus reasonable to treat as an analogue of 
the entanglement encountered in systems of distinguishable particles only an additional amount 
of correlation going beyond that stemming from symmetry requirements. 

We proposed a way of treating all non-classical correlation, i.e., those which can be identified 
with the 'genuine entanglement' and not caused merely by symmetries. This unifies all cases: 
of distinguishable particles, fermions, and bosons, and can be easily extended to hypothetical 
multipartite systems consisting of particles subjected to arbitrary parastatistics. 

We defined simple (non-entangled) pure states as one-dimensional selfadjoint projectors as- 
sociated with simple tensors obeying appropriate symmetries identifying particles as bosons, 
fermions, etc. Consequently, simple (non-entangled) mixed states were defined as convex com- 
binations of simple (non-entangled) pure ones. Such an unifying approach allowed also for 
description of such tools, known from the entanglement theory, as the Jamiolkowski isomor- 
phism and Schmidt rank, to systems with other symmetries. The introduced concept of S-rank 
not only provides us with a tool for distinguishing entanglement of pure states with a given 
parastatistics, but is interesting also per se, as it offers the simplest characterization of highest 
weight vectors we know. 

In the case of two fermionic subsystems our approach identifies non-entangled pure states to be 
the same as in all other approaches mentioned in Introduction, i.e., we identify them with simple 
antisymmetric tensors in the meaning explained in Section [3] above. In the bosonic case, from 
the geometric point of view, we clearly have a priori two inequivalent types of non-entanglement: 
tensor products of identical states and states associated with symmetrizations of products of 
orthogonal vectors. Non-entangled states of two different types are not connected by local 
unitary transformations which is in contrast to the familiar situation of distinguishable particles 
and intuitions build upon the fact that all separable states of distinguishable particles can be 
obtained from a single one by local transformations. Although this is obviously acceptable, it 
poses an open fundamental problem what is a physical meaning of two geometrically inequivalent 
types of non-entanglement. 

In our approach we adopted the view that non-entangled pure bosonic states are simple sym- 
metric tensors - tensor products of identical vectors. We find at least two arguments justifying 
this choice. In |2] it was pointed that all states which are symmetrizations of products of dis- 
tinct vectors can be used to perform such clearly 'non-classical' tasks like teleportation. This 
definitely remains in conflict with the basic intuition connecting non-entanglement with the 
purely classical world. Second, from purely mathematical point of view, only tensor product 
of identical vectors provide highest weight vectors of the corresponding representation of the 
unitary group, like in the other cases. 

Another achievements of our paper are explicit characterizations of simplicity (non-entanglement) 
for pure states in terms of directly verifiable, quadratic in coefficients, conditions which are com- 
putationally much easier than those proposed in the literature. 
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